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Abstract. We study zero-temperature XX chains and transverse Ising chains and 
join an initially separate finite piece on one or on both sides to an infinite remainder. 
In both critical and non-critical systems we find a typical increase of the entanglement 
entropy after the quench, followed by a slow decay towards the value of the 
homogeneous chain. In the critical case, the predictions of conformal field theory 
are verified for the first phase of the evolution, while at late times a step structure can 
be observed. 



1. Introduction 

The complexity of many-body quantum states manifests itself in the connection between 
different parts of a system in the wave function. It can be measured by the entanglement 
entropy S which has been studied extensively over the last years and shows, for example, 
a logarithmic divergence at quantum phase transitions in chains [I]. More recently, its 
temporal evolution for evolving quantum states has come into the focus of the studies. 
The simplest situation is a quench, where the Hamiltonian is changed instantaneously 
from Hq to Hi and one follows the subsequent evolution of the initial state, usually 
taken to be the ground state of Hq. The first such studies considered global quenches, 
in which a parameter is changed everywhere in the same way [21 13] . It was found that 
in this case the entanglement entropy becomes extensive, in contrast to the equilibrium 
states where it is proportional to the interface area between the two chosen parts of the 
system. 

A different situation is that of a local quench, where a parameter is changed only in 
one or in a few places. For example, a single nearest-neighbour bond could be modified. 
In one dimension, this can have strong effects, since by removing or adding the full 
bond, one cuts the chain or joins previously separate pieces. In the first case, the 
initial entanglement is conserved, but in the second case it shows interesting temporal 
behaviour. This was investigated recently in a numerical study of free electrons hopping 
on a chain with initially one defect in form of a weakened bond which was then removed 
[4]. In spin language, this corresponds to an XX model. The entanglement entropy 
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was found to increase to a maximum, or plateau, followed by a slow decay to its 
equilibrium value. The maximum value was proportional to InL, where L is the size 
of the chosen subsystem, and the whole plateau could be described by a logarithmic 
expression containing L and the elapsed time t. In a subsequent publication, Calabrese 
and Cardy [5] showed how to obtain such plateaus and the corresponding formulae from 
conformal field theory. They also found that the long-time tails in S are absent in this 
continuum treatment and must therefore be lattice effects. 

In the present study we continue these investigations by looking at chains where a 
finite segment is initially fully separated from the rest and then joined to it. This can be 
done on one or on both sides of the segment and the rest is assumed to be infinite. All 
parts of the total system are initially in their ground states. As examples we consider 
the XX model and the transverse Ising (TI) model for which the entanglement can be 
obtained in a relatively simple way from the correlation functions. This is outlined in the 
following section 2. The critical case is then studied in section 3. As for the single defect, 
we find logarithmic plateaus for S{t) which turn out to be in very good agreement with 
conformal predictions. For the doubly-infinite geometry we derive the corresponding 
formula by extending the treatment in [5]. We also illustrate how the entanglement 
process manifests itself via fronts in the eigenvectors of the reduced density matrix 
and show the evolution of the single-particle spectra underlying S. An interesting new 
phenomenon appears in the long-time behaviour of S{t), where one finds a step structure 
in the slow decay towards equilibrium. The non-critical case is investigated for the TI 
model in Section 4. Here we find a plateau in S as well, but with a flat shape (up to 
oscillations) and a limited height depending on the transverse field. Our findings are 
summarized in Section 5. Two appendices treat the relation between the entanglement 
entropies of the XX and the critical TI model and the derivation of the conformal 
formula for one of the geometries. 



2. Models and geometries 



In the following we study the XY quantum chain with Hamiltonian 
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where the a" are Pauli matrices, for two special cases, namely 

• K, = h = which is the XX chain, critical and equivalent to a half-filled electronic 
hopping model 

• K = 1 which is the Ising chain in a transverse field. 

As is well known, the operator H becomes a quadratic form when expressed in 
terms of Fermi operators c„, cj^ and reads after diagonalization 

^ = E^.(*.- 1/2) (2) 
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with 

Uq = \J sin^ q+ {h + cos g)^ . (3) 

In general, the rjq are obtained from the c^, via a Bogohubov transformation with 
two sets 4>q{n), ipq{n) of coefficients [6]. For the XX case this reduces to a simple linear 
combination of the c„ and one can also write H as 

^xx = -^cosgr/Jr/g (4) 

with q = Tm/{N + l),n=l,2,...A^ for open boundary conditions. 

We will investigate the two situations depicted in Fig. [H In the first case the system 
consists initially of a segment of length L and a rest of length N — L where N ^ L, 
both of which are in their ground state. We will call this the semi-infinite geometry. In 
the second case, the segment is sandwiched between two such environments. This will 
be called the infinite geometry. At time t = 0, the missing bond(s) will be added and 
the time evolution of the state monitored. 



Figure 1. The two geometries studied in the paper. 

We will be interested in the entanglement of the segment with the rest of the system. 
The corresponding entanglement entropy 5* follows from the reduced density matrix p 
of the segment which has the form [7j 

P = ^e-«, n = jZeu{t)flh (5) 

k=l 

Here Z is a normalization constant ensuring Tr p = 1 and the fermionic operators fk 
follow from the Cn by an orthogonal transformation. Then S = — Tr (p In p) is determined 
by the single-particle eigenvalues Skit) according to 

S{t) = - E (kit) \nCk{t) - E(l - Ck{t)) HI - Ck{t)), (6) 

k k 

where Cfc(t) = l/{exp{ek{t)) + 1). 

Since the initial state is a Slater determinant, the ek{t) resp. the Ck{t) follow from 
the one-particle correlation functions at time t. These can be combined into a 2N x 2N 
matrix 

/ (F^ F^^ (F^ F^^ \ 
\ (F2 Fi) (F2 F2) j 

where F^ = (cj^ + Cn) and F^ = — ^(c.|^ — c„). Restricting this matrix to the L sites in the 
subsystem, its eigenvalues are Ca; and 1 — Cfc- In the XX case, it is sufficient to consider 
the N X N correlation matrix 

C = (4c„) (8) 
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which has eigenvalues (k when restricted to the subsystem. 

The time evolution of the fermionic operators with the final Hamiltonian gives the 
correlation matrix at time t in terms of its initial value as [8], [9] 

G(t) = R(t)G(0)R^(t) (9) 

where R(t) is a rotation matrix with elements which are sums containing the (l)q{n),iljg{n) 
and exponential factors exp{iujqt). In the critical case and taking the thermodynamic 
limit, they can be expressed in terms of Bessel functions. The initial correlations are 
particularly simple for the XX model, where one finds for a segment with sites between 
1 and L 

sin(f(2-j)) _ sin(f(^ + j)) 1 
.sin(2(fTT)(^-^')) sin(^(z + j))J 
It was shown recently that in the ground state there is a close connection between 
the entanglement entropies of the XY chain with h = and the TI chain [10]. The 
consideration can be extended to the time-dependent case as sketched in Appendix A 
and gives at criticality 

Sxx{L,t) = 2STi{L/2,t/2) (11) 

One can see this relation clearly in numerical computations of S. For the spectra it 
means that the same Ek occur in both systems, which one can also verify numerically. 
In the critical case we will therefore only present results for the XX model. These 
were obtained by determining Cij(t) in the subsystem by summing Bessel functions as 
in [Ij. For the non-critical case we will consider the TI model. Here the procedure 
was to calculate the full correlation matrix G{t) from G(0) for a sufficiently large total 
system using ([9]) and then to specialize to the subsystem. The size was chosen such that 
refiection phenomena in the larger subsystem could not yet enter the dynamics. 



3. Critical chains 

General behaviour 

The general behaviour of the entanglement entropy is shown in Fig. [2] and can be 
described as follows. Starting from the value 5(0) = (since the subsystems are 
separated initially), it rises quickly to a broad plateau which extends to the times 
t = 2L in the semi-infinite and t = L in the infinite geometry, respectively. Afterwards it 
decreases slowly to its equilibrium value. The characteristic times are those an excitation 
with the maximum velocity vp = ^ needs in order to travel from the initial gap to the 
end of the segment. In the semi-infinite case this involves a refiection at the boundary, 
which gives t = 2L. 

Such a propagation phenomenon is clearly visible in the eigenvectors of the single- 
particle states of p. Fig. [3] shows the lowest one for four different times and both 
geometries. In the semi-infinite case (left side) one sees a front propagating from the 
initial defect towards the end where it is refiected. For the infinite geometry (right 
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infinite 




Figure 2. Entanglement entropy for L = 40 and both geometries. The dotted hnes 
on the right are the corresponding equihbrium values. 



side) one has two fronts propagating inwards, penetrating each other and continuing 
afterwards until they reach the ends of the subsystem. 




10 20 30 40 50 60 70 80 90 100 10 20 30 40 50 60 70 80 90 100 



Figure 3. Front propagation in the lowest single-particle eigenvector for L = 100, with 
time increasing from top to bottom. Left: semi-infinite geometry, t = 20, 60, 120, 160. 
Right: infinite geometry, t = 10,30,50,70. Shown are the absolute squares of the 
amplitudes. 

The evolution of the single-particle eigenvalues is shown in Fig. |H The horizontal 
scale has been chosen in order to make the features at both small and large times clearly 
visible. The general picture is similar to the one that was found previously for one 
defect in an infinite system [4j. In the beginning the Ek drop to some transient values 
which produce the plateau in S{t) and subsequently they approach their equilibrium 
values. While this happens fast for the lowest Ek, considerable deviations remain for long 
times among the higher ones. These are the "anomalous" eigenvalues. The deviations, 
however, are stronger than in and there is a new feature in the form of certain 
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intervals where the eigenvalue is roughly constant before evolving again. We will come 
back to this phenomenon below. 
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Figure 4. Time evolution of tlie lowest four single-particle eigenvalues for the semi- 
infinite geometry with L — 40. The parameter to was chosen such that t = 2L gives 1 
on the horizontal axis. 



Plateau region 



For the semi-infinite geometry, Calabrese and Cardy recently derived the following 
formula for S{t) from conformal field theory [5] 



S{t) 



6 



In 



4L 



-t sin 



L TT 



7rt\ 
2l) 



+ k'. 



(12) 



Here c is the central charge, equal to 1 for the XX model, and the cutoff parameter 
e appearing in the calculation has been included in the constant k'. This formula is 
supposed to be valid as long as the sine remains positive, i.e. for times up to t = 2L. 
The lower curves in Fig. [5] show a comparison of numerical data with this prediction for 
a system wit L = 60. There are deviations at the right and left end of the interval and 
the data shows additional oscillations, but overall the agreement is very good. Leaving 
the value of c free, a fit between t = 10 and t = 110 gives c = 0.993. Also the dependence 
on L is described properly, i.e. curves for different L collapse in reduced variables. Thus, 
in spite of an approximation involved deriving it [5j, formula (fT2!) describes the lattice 
result very well. 

The infinite geometry was not treated in [5]. However, as sketched in Appendix B, 
one can generalize the approach to this case by using an additional conformal mapping. 
The result after a rather long calculation is 



Sit) 



In 



2L 



-t sin 



TT 



TTt 

T 



+ k. 



(13) 



Thus, compared to (fT2l) . the prefactor has twice its value, which is due to the two points 
of contact between subsystem and surrounding, and 2L has been replaced by L. A 
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Figure 5. Entropy plateau. Numerical results (solid) and conformal predictions 
(dashed) for L = 60. Upper curves: infinite case, t — t/L. Lower curves: semi-infinite 
case, T = t/2L. 



comparison of this formula to the numerical data is given by the upper curves in Fig. 
ini again for L = 60. The curves look very similar to the lower ones and the agreement 
is again very good. The same kind of fit as before gives in this case c = 1.007 

The conformal field theory, which deals with particles of uniform velocity v = 1, 
does not give the oscillations seen in the numerics and the relaxation at long times. 
Both are lattice effects connected with the dispersion of the elementary excitations in 
the model. 

Long times 

The behaviour at long times is more intriguing here than for a single defect in an infinite 
chain. From Fig. [2] it might seem that S{t) simply decays towards its equilibrium value. 
Looking closer at the semi-infinite case, however, one finds a clear step structure. This is 
shown in Fig. [6] for two values of L. These steps become quite regular if one plots 5* vs. 
1/t as in the right part of the figure. Their height decreases with L roughly according 
to a power law L^" where a ~ 0.90. A comparison with the single-particle spectrum in 
Fig. m shows that the steps occur for those times where the anomalous eigenvalue has 
the fiat intervals mentioned previously. These are short for the earlier times and longer 
for the later times which gives the regular structure in 1/t. 

If one examines the eigenvector of the anomalous single-particle state one finds 
another surprise. For the fiat intervals, it becomes remarkably simple. This is shown 
in Fig. [7| for L = 60 and two times related to the third-last and the last plateau, 
respectively. Plotted are the absolute squares of the amplitudes. These pictures suggest 
the functional form 

\(f){j)\^ = Asm\qj) ; Q = J^^^ (14) 

as for a standing wave in the isolated subsystem and the comparison shows that this 
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Figure 6. Entanglement entropy at large times for the semi-infinite geometry. Left: 
5 vs. t. Right: S* vs. 1/t. The equilibrium value Sq is subtracted. 



indeed fits the data very well. Moreover, one can take these g-values and calculate the 
time an excitation with Uq = — cos q created at the initial defect needs to travel through 
the system. With Vg = sin q this is 

T^^L^LiL±}l (15) 

Vn TT n 

The times so calculated fit the plateau times in ek{t) resp. S{t) remarkably well. The 
expression also explains the scale used in the right part of Fig. [61 The plateaus then 
appear at integer values of the variable. 
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Figure 7. Anomalous single-particle eigenvector for L = 60. Left: t = 390. Right: 
t = 1160. The dotted lines are fits with the function 



Taking all together, one can say that the discrete levels of the isolated subsystem 
reappear in the time evolution. Since for long times the smallest g-values enter, one is 
seeing here (in contrast to the fronts discussed earlier) the slowest excitations in the 
subsystem. However, the selection of the g's is not by interference of a particle with 
itself, since (1151) gives only the time for travelling one way. 

All considerations so far were for the semi- infinite case. In the infinite geometry, 
one also finds such plateaus in the Ekit) at times given by (ITSll with L replaced with 
L/2. The corresponding eigenvector also has the same features. However, one does 
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not observe plateaus in S{t), because they are masked by the contribution of another 
non-constant eigenvalue nearby. 

4. Non-critical chains 




I ' ^ ' ^ ' ^ ' ^ ' 1 

50 100 150 200 

t 

Figure 8. Entanglement entropy for the TI chain in the semi-infinite case with h = 2 
and L = 20. The dashed horizontal lines indicate the plateau value and the equilibrium 
value, respectively. The arrows give the four longest times r„ = L/vn defined in the 
main text. 

We now present results for the non-critical TI model in the disordered region h > 1. 
The general behaviour of S{t) is shown in Fig. [8] for the semi-infinite geometry and 
L = 20. This is a relatively small size but still much larger than the correlation length 
^ = {h — 1)^^ for the chosen value of h. As before, one sees an initial plateau which 
extends up to t = 2L. This is consistent with the velocity of the elementary excitations 

Vq = h sin q/uq (16) 

which gives a maximum value v = 1 for g = 7r/2 as in the critical case. The height of the 
plateau is independent of L for L ^ ^ and decreases monotonically as one moves away 
from the critical point. At the shown value /i = 2 it is already rather small. Its detailed 
behaviour can be seen in Fig. [91 One finds numerically that the plateau height is very 
accurately given by ^/^TSo where So{h) is the equilibrium value of the entanglement 
entropy for large L given analytically in [11]. The relative deviation from this formula 
is less than 0.1%. 

On top of the plateau one sees oscillations as in the critical case. Due to the 
small plateau height, however, they are more prominent here. They reflect the lattice 
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Figure 9. Plateau values of the entanglement entropy as a function of the field h. 
The line is the equilibrium value So{h) multiplied by ^/^T. 



structure of the chain and are shown in more detail on the left side of Fig. [101 Their 
amplitude decays as 1/t while their period changes with the transverse field. From the 
Fourier spectrum on the right side of the figure one sees that they contain two main 
frequencies. These are given by u!± = {h + 1) ± {h — 1) where h + 1 and h — 1 are 
respectively the upper and lower band edge of the dispersion curve ([3]). They would 
appear in a stationary phase approximation for the integrals in the rotation matrices 
R(t) for the time evolution. Of these, the larger frequency a;+ = 2h contributes with 
higher weight. Moreover, there is a peak in the spectrum at zero frequency. Thus, to 
leading order S(t) in the plateau region is described by the expression 



S{t) ~ V^Soih) + + Bcos{2ht)] 
where A and B are constants depending a priori on L and h. 



(17) 
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Figure 10. Left: Entanglement entropy in the plateau region for L = 30, h = 1.4 
(top) and h ~ 2.6 (bottom). Right: Power spectrum of t[S{t) — y^So] for h — 2.6. 
The main peak is located at to ~ 5.2. 



For times longer than 2L, Fig. [8] shows a slow decay of the entanglement entropy 
which is similar to the one in the critical case and the curve seems to approach the 
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equilibrium value Sq. If one assumes that, the decay is approximately a power law 
with an exponent near one. Due to the calculational procedure, the accessible times are 
smaller here and limited to about 300. In this region, the spectrum of the single-particle 
eigenvalues still deviates markedly from its asymptotic form. Only the lowest one has 
the value one finds from the relation to corner transfer matrices [T2l ITT] , namely 



where I{k) is the complete elliptic integral of the first kind, while the next two are 
irregular. Thus one cannot check the approach to equilibrium very well. The curve for 
S(t) shows again certain structures, but these are nodes between oscillatory regions 
rather than steps. Nevertheless, as can be seen from Fig. [HI they occur at the 
characteristic times T„ = L/vn calculated from the velocities (|T6l) . The allowed values 
of q in the non-critical case are not equidistant and lead to different velocities near the 
upper and lower band edge. Finally, one should note that the overall behaviour of S{t) 
has a certain universal character, since only one time scale, namely L, enters. Thus one 
obtains an approximate collapse of the curves for different L (but fixed h) if one plots 
them as a function of t/L. This means that after times t > 2L, the decay is slower for 
larger subsystems. 

5. Summary and discussion 

We have studied a special case of a local quench, namely the joining of initially separated 
parts of a one-dimensional system, all prepared in their ground states. The quantity of 
interest was the entanglement entropy between a finite piece and the remainder of the 
chain. We found that S behaves in a similar way as after the removal of a single defect 
from an infinite system [T]. There is a transient regime where entanglement builds up 
rapidly and a plateau forms. Its length is determined by the velocity of the fastest 
particles. In the critical case, its form is given by a universal logarithmic function and 
follows from conformal field theory. For non-critical chains, it is rather fiat and its 
height varies with the parameters of the system. At later times, there is a slow decay 
towards the equilibrium entanglement in all cases. Thus one has a kind of overshooting 
phenomenon. 

The critical case turned out to be particularly interesting for the semi-infinite 
geometry. Then the final decay of S proceeds in a step-like fashion with fiat regions 
determined by the travelling times of the slowest particles in the initially isolated 
subsystem. Thus, while the gradual descent is generally a consequence of the dispersion 
in the velocities, it shows here some memory of their initial discreteness. This is 
connected with a special behaviour of the anomalous single-particle eigenvalues of the 
reduced density matrix and their eigenvectors. On that level, it seems rather general 
for this type of quench. However, to be observable in the entanglement entropy, a 
proper eigenvalue structure is necessary. Although we can describe the phenomenon 
rather precisely, we have no complete picture of the mechanism. The same holds for the 



ei = 7cl{k')/l{k), k = l/h, k' = VT^^ 
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simple formula giving the plateau height for the TI chain in terms of the equilibrium 
entanglement. 

The physical situation we have addressed, has been considered with a different 
focus in various other papers. For example, in [13] oppositely magnetized XX half- 
chains were brought into contact and the interest was in the magnetization profile. In a 
fermionic picture this corresponds to charge transport via tunnel contacts as studied in 
[H] . Generalizations include the case of finite temperatures and a sandwich structure 
corresponding to our infinite geometry [151 HBl 113, EH] • In [19] also the entanglement for 
large times was considered. However, for finite temperatures the entanglement entropy 
becomes extensive and does not measure the quantum properties we were studying here. 
Nevertheless, such an extension can put the zero-temperature features in a broader 
context. 
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Appendix A: Entanglement in XY and TI chains 

As pointed out recently [10], there is a close connection between the ground state 
entanglement entropy of the zero field XY chain and the TI chain. This was derived 
by a detailed comparison at the level of the fermionic Hamiltonians and the correlation 
matrices. Here we extend this to the non-equilibrium case using more general arguments. 

It has been known for a long time [201 EI] that the zero field XY Hamiltonian ([1]) 
can be split into a sum of two independent TI Hamiltonians by distinguishing even and 
odd lattice sites and using a dual transformation. If the couplings (1 ± k) in ([T]) are 
denoted by and Jy these are 

N/2-1 ^ N/2 

Hi, I = - 9 E J^n'^n+i - 9 E Jy^'^ (19) 
^ i=i ^ i=i 



and 



N/2-1 N/2 
-L - x,2 x,2 J- v-^ 7 2,2 



Hi,2 = -2 E Jy^' <^ - 2 E J^i (20) 

1=1 i=l 

Here the r operators are related to the a's via [lOj 

2i-l 2j-l 



j=i i=i 
^^'1 — ^y -^^'2 



2,1 y y z,z X X /oi \ 

= f^2i-l4i , = ^2i-l(^2i (21) 

One then has the decomposition 

HxY = I [Hi,i + Hj,2] (22) 
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with [Hj^i, Hi^2] = 0. If the XY chain is cut at one point, the TI chains are also cut. 
Therefore, the initial ground state of the XY chain factorizes as 

$xy(0) = $/,i(0) $7,2(0) = <l>},i(0) 1(0) $1^2(0) $?,2(0) • (23) 

where the superscripts refer to the two subsystems. The evolution operator Uxyit) = 
Q-itHxY factorizes as well, since [Hj^i,Hj^2\ = 

Uxrit) = Ui,i{t/2) f//,2(t/2) (24) 

and the state at a later time is given by 

<l>xy(t) = $/,i(t/2)$7,2(t/2). (25) 

Tracing over the unwanted degrees of freedom, one obtains the reduced density matrix 

as 

Pxvit) = Pi,i{t/2) p/,2(t/2) (26) 

and, restoring the length L, one has for the entanglement entropy 

Sxy{L, t) = 5,,i(L/2, t/2) + Sj^2iL/2, t/2) . (27) 

In the isotropic case, i.e. for the XX chain, the two TI models are identical and critical 
and (ET]) gives in Section 2. 

Appendix B: Conformal mapping 

For the semi-infinite chain, the formula f|T2l) was derived in [5] within a path-integral 
approach for the quantum system. This leads to considering a complex (x, ir) plane, 
where r is the imaginary time. In this plane, one has a vertical line representing the 
boundary and a pair of half-infinite vertical lines representing the cut in the chain which 
exists up to ±ie. The entanglement entropy S{t) is obtained from the expectation value 
of an operator $„(zro) in the gap between these half-lines after the analytic continuation 
To — >■ it. To obtain this expectation value, one first simplifies the geometry by mapping 
the area to the right of the boundary line, where the path integral is to be taken, onto the 
right half-plane. As explained in [5], this can be done approximately by the combination 
of a Joukowsky transformation, a reflection at the unit circle and a logarithmic map. 
The resulting transformation is analytically invertible which is needed for the application 
of the conformal formulae. 

The infinite chain with initially two cuts can be related to this problem in the 
following way. One has to consider in this case a complex plane with two pairs of half- 
infinite vertical lines and the path integral is taken everywhere. The situation is shown 
in Fig. [11] on the left. The transformation 

z' = ^iVz^ + - z) - L) (28) 

which is a Joukowsky transformation as in [5], followed by a reflection, an inversion 
and a translation by L then maps this approximately onto the shaded region of the 
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Figure 11. Conformal mapping of tlic infinite geometry with two cuts (left) to the 
semi- infinite one with one cut (center) and onto the right half-plane (right). The path 
integrals are taken in the dark areas. The colours show how the various lines transform. 



2;'-plane shown in the center of Fig. [TTl This is exactly the geometry for the single 
cut and one can use the previous procedure to map it to the right w' half-plane. The 
only additional feature is that the entanglement entropy follows here from a two-point 
function ($„(zro)$„(z'ro — L)). From its form in the w'-plane [5] 



($„(^)$„(w^)) 



(29) 



\w[ - \2Rew[\ \2Rew'2\ 

where Xn = c{n — l/n)/12 one obtains it in the z-plane by the usual conformal 
transformation formula. 



dw' 



dz 



dw' 



dz 



2 J 



(30) 



After a lengthy calculation and analytical continuation, assuming t,L ^ e, this leads 
to 

.2Lj ' " ' ^ 



(31) 



J J sin^ (Tct/L) 

Here the first two factors in the brackets come from the derivatives dw' /dz while the 
sine arises from the denominator in (|29l) . The entanglement entropy is the negative 
derivative with respect to n at = 1 which gives f|T3|) in section 3. 
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